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Abstract 



The coupling of the baryonic current to new neutrino states v^, with strength in excess of 
the weak interactions is a viable extension of the Standard Model. We analyze the signature 

{Sj ' of z/f, appearance in the solar neutrino flux that gives rise to an elastic scattering signal in 

dark matter direct detection and in solar neutrino experiments. This paper lays out an in- 

^ ■ depth study of i/& detection prospects for current and future underground rare event searches. 

We scrutinize the model as a possible explanation for the reported anomalies from DAMA, 
CoGeNT, and CRESST-II and confront it with constraints from other null experiments. 



1 Introduction 

The phenomenon of neutrino mass-mixing and thereby induced lepton-flavor oscillations 
constitute the only conclusively detected deviation from the Standard Model (SM) to date. 
Though the neutrino sector is certainly the most elusive part of the SM, a tremendous 
and sometimes painstaking experimental effort has firmly established the patterns of mass 
splitting and mixing for the three SM neutrinos, u e , z/ M , v T . The multitude of data is con- 
sistent with neutrino flavor eigenstates being the linear combination of — at least — three 
massive ones v\, u 2 , 1/3, with mass-squared differences Am^i = (9(10 _5 )eV 2 and Am^j ~ 
Am| 2 = ±(9(10~ 3 ) eV 2 and mixing angles with magnitudes sin 2 26 , 13 ~ 0.1, #23 ~ 7r/4, 
tan 2 8 12 ~ 1/2 [IJ|2]. Neutrinos are predominantly detected via their charged current (CC) 
interaction on matter, and elastic scattering (ES) on electrons. The only pure neutral cur- 
rent (NC) process for solar neutrinos with participation of all active neutrino flavors observed 
remains the deuteron breakup reaction d + v x — > p + n + v x in the SNO experiment [H H] . 

On the theoretical side, the interest to the elastic scattering of neutrinos on nuclei dates 
back to Drukier and Stodolsky [5], who outlined their vision of a true neutrino observatory: 
neutrinos with MeV-scale energies as they emerge from the interior of the sun, from nu- 
clear reactors, from spallation beam experiments, or from supernovae explosions can scatter 
elastically and with N 2 enhancement in the cross section on keV-scale recoiling nuclei of 
N neutrons [6]. Even though the energy transfer to the target nucleus is small, the idea 
was that its effect in a detector at cryogenic temperatures may be macroscopic given that 
the specific heat of the material can be minute. The underlying principle of NC neutrino 
scattering on nuclei has found its proliferation in Dark Matter (DM) searches for weakly 
interacting massive particles (WIMPs) [7j. 

Direct detection experiments seek for evidence of DM via its elastic scattering on various 
targets such as crystals made from germanium or liquefied nobles gases such as xenon. No 
univocal evidence for DM has yet been found but upper limits as low as 7 x 10~ 45 cm 2 [S] in 
the WIMP-nucleon cross section have been inferred for the optimal range of recoil energies. 
These experiments have reached a level in sensitivity such that neutrino coherent scattering 
on nuclei is being discussed as potential background for future ton-scale experiments; see 
e.g. HD|. The latter process has a cross section a u /N 2 ~4x 10~ 43 (£' I/ /10MeV) 2 cm 2 
where E u is the neutrino energy. Indeed, the average flux of neutrinos at earth, dominated 
by the solar pp-chain of thermonuclear reactions, $ pp ~ 6 x 10 10 cm _2 s _1 [ll\, exceeds 
the flux expected from WIMPs, $dm ~ 10 5 (100 GeV/mDivi) cm" 2 s _1 , by many orders of 
magnitude. The reason why neutrinos are not copiously detected via their NC interaction 
in DM experiments lies in the soft recoil spectrum they induce, E^ < 2E 2 /m^ ~ 2keV x 
(100/A) (E u /10MeV) , where m^ and A are the mass and atomic mass number of the 
target nucleus. Whereas these experiments fall short in sensitivity to SM neutrinos, they 
may nonetheless be powerful probes of an extended neutrino sector [12J. 

This paper surveys potential signatures of the existence of yet another "sterile-active" 
neutrino state, z/f,, where active is understood in the sense that v^ shall interact via a new 
neutral current interaction with baryons (NCB) [12]. This "baryonic" neutrino is also "ster- 
ile" in that it shares no SM NC or CC interactions. In particular, this implies that vj, does 
not scatter elastically on electrons. Remarkably, v h can then be coupled to baryons with a 
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strength Gb which exceeds that of the weak interactions substantially, Gb/Gf = (10 2 — 10 3 ); 
Gp is the Fermi constant. This, in turn, implies the existence of a force mediator with the 
mass much smaller than mw The key observation made in [T2] is that for MeV-scale energies 
of Ub the ratio of elastic to inelastic cross section with nuclei scales as 

^(elastic) ^ 8 (A_\ 2 /10MeV\ 4 /2_fm\ 4 
a VbN (inelastic) ~ \iw) \ E v ) \ R N J ' U 

where Rn denotes a nuclear radius. This tremendous enhancement opens an exciting phe- 
nomenological window of opportunity: v\, states can be searched for in the growing number 
of direct detection dark matter experiments. The less constrained choice of parameters for 
this model is when the oscillation length for Vqm — > v b transition is long enough not to create 
any significant fluxes of reactor/beam/atmospheric neutrinos, while for solar neutrinos one 
can have a sizable fraction of v\>- In particular, solar 8 B neutrinos have the best combination 
of large flux and high end-point energy for producing an observable signal in rare event 
searches. Relation ([I]) makes even small scale DM searches competitive with dedicated large 
target-mass neutrino experiments in their sensitivity to v^ that originate from 8 B neutrinos. 

Moreover, as proposed in [12], if the oscillation length of Uj, is on the order of the earth- 
sun distance, the signal may be annually modulated in a non-trivial way, with a possibility 
of a bigger flux in the summer, and with modulation amplitude larger than naively expected. 
In this paper we analyze this corner of the parameter space in great detail, as it could offer 
a //^-scattering explanation to the DAMA signal [131 03] while still being compatible with 
other null-results from DM experiments; for a relevant collider study in this context cf. [15]. 
The long-standing DAMA "anomaly" has very recently received some additional impetus by 
the reports [TH], [TT| of the CoGeNT collaboration on 1) an unexpected rise in observed events 
at low nuclear recoils and 2) an indication of annual modulation, see e.g. [T7J [T5J [THJ, |2"U] . 
Finally, the CRESST-II experiment has published results of its latest run, which had some 
population of unexpected events on top of existing backgrounds [21] . 

The attempts to explain positive signals (DAMA) and possible hints on non-zero sig- 
nals (CoGeNT, CRESST) while staying consistent with null results of other groups, are 
widespread in the WIMP literature, see e.g [221 1231 122 ESI M, 123 123 1221 EQ] and references 
therein. The models that fare best feature ~ 10 GeV WIMP masses, although the overall 
consistency of the "light WIMP" explanation for the "DM anomalies" remains doubtful. In 
this paper, we provide an in-depth critical assessment of v^ models with regard to their poten- 
tial to explain positive signals, and make predictions for the upcoming experiments, paying 
attention to those that could potentially differentiate between v\, and light dark matter re- 
coils. We leave other phenomenological aspects of this interesting, but to date a relatively 
poorly explored model to subsequent work. 

This paper is organized as follows: in the next section, the baryonic neutrino model is 
reviewed. In Sec. [3] we cover the current and future sensitivity of DM searches to i/& and 
study the potential explanation of the DAMA, CoGeNT, and CRESST-II signals. In Sec. H] 
the elastic scattering of v\, in neutrino experiments is considered, and in Sec. [5] we reach our 
conclusions. 



2 Baryonic neutrinos 

In the baryonic neutrino model [12] , the SM gauge group is extended by an abelian factor 
U(l)s. The new neutrino is a left-chiral field v b = |(1 — 7 5 )z/& with charge q b = ±1 and 
gauge coupling gi > 0. Leptons are neutral under U(1)b but all quarks q = Ql, Ur, d,R carry 
baryonic charge 1/3 with gauge coupling g b > 0. The SM Lagrangian is supplemented by 
the following terms 

C B = vrfiidp - giq b V,)u b - l - gb J2 cfW - \v^v^ + \m 2 v V^ + C m . (2) 

g 

We have assumed that the new gauge field V^ with field strength V^ u has acquired a mass m v 
by the spontaneous breakdown of the U(l)# symmetry with Higgs^ VEV {(p b ) = V{,/y/2; the 
sum in the second term runs over all SM quarks q. The part C m is responsible for generating 
neutrino masses and the mixing between SM neutrinos and the new state v b . A simple 
UV-completion of C m is one where — once electroweak symmetry and U(l);, are broken — new 
right-handed neutrinos vr induce mass terms for SM neutrinos as well as for v b . 

(v' L \ /0 ml 0\ 

C m = -Nl&MN L + h.c, N L =\vg\, M=\m D m R v b b . (3) 

2 V'J V° v b b T Oj 

Here, v' L denote the three SM neutrinos with Dirac mass matrix mo, v'^ = Cv R are 
the charge conjugate states of v' R with Majorana mass matrix rriR, and b is a vector of 
Yukawa couplings generating mass for v' h . In the simplest case the new right-handed neutrinos 
generate masses for u' L and v' h simultaneously. Introduction of a right-handed partner for 
v' h will cancel a U(l)^ triangle anomaly and remaining gauge anomalies can be cured by 
the introduction of a new family of heavy fermions with appropriate quantum numbers [31J. 
From the phenomenological viewpoint the details of this will not be of importance for this 
work. 

The mass matrix M in ()3]) is diagonalized by A^diag = {Yl) MV£ where V£ is a unitary 
matrix, defining the mass eigenstates Nl = {V^)'N' L . Diagonalization of the charged lepton 
mass matrix by 3 x 3 unitary matrices V[ R with mass eigenstates II = (V[] )H' Ll and Ir = 
(V R yl R where I = (e~,fi~,r~) then determines the mixing among the SM active neutrino 
states. Assuming that the eigenvalues of rriR are much larger than any other values in M, 
the seesaw mechanism is operative and one can integrate out the heavy, right handed states. 
We are left with the 4x4 mixing matrix U which connects "flavor" v a i (a = e, /i, r, b) and 
mass rikL (k = 1, • • • , 4) eigenstates: 



e 



12 3 4 

( -\ 
b V • • -j 



(4) 



where Upuns is the usual 3x3 mixing matrix among the active flavors pp. The NCB current 
in the respective interaction and mass eigenbasis reads 

3 ncb = v bl^Vb = ^2 u lkU±k>n k Ll^n k , L . (5) 

k,k' 

2.1 Neutrino oscillations and matter effects 

Aiming at a scenario in which the baryonic neutrino is coupled stronger than Gp, the ques- 
tion which immediately arises is whether new matter effects are to play a role. The index of 
refraction of ^-propagation in matter is found by computing the forward scattering ampli- 
tude of u b on matter, described by the effective Lagrangian 

£ cff = -GbJncb E W ^ N ' Gb = «• ^T = ^°f- (6) 



N=n,p 



m 



v 



In the last equality we have introduced the parameter M > to measure Gb in units of 
Gp with the sign of the interaction determined by the charge q b of v b . In an unpolarized 
medium, one obtains the following matter potential 

V NC B = ±q b MG F n B (Y N + 2Y Ub ), Y f = -1 f -, (7) 

tib 

where the plus sign is for v\> and the minus sign for X7 & ; Yf is the particle-antiparticle asym- 
metry, normalized to the number density of baryons Ub- The first term in the first equation 

describes the potential in ordinary matter while the second term is the potential for v hL in 
a hypothetical sea of baryonic neutrinos. 

In ordinary matter, with mass fraction X p in form of bound or unbound protons, the 
induced matter potentials (up to coherence factors) compare as follows 

V NCB : Vcc : V NC = QbM : V2X P : -V2(l - X p )/2, (8) 

where we have made use of the charge neutrality condition. As is evident, for M ^> 1 
baryonic neutrinos experience the largest matter effect since X p is always of order unity. 

The large strength of the NCB interaction may suggest that the flavor evolution in matter 
is dominated by Vncb- In a simplified two- neutrino case the Schrodinger-like equation 
describing the transition probabilities P a ^/3(x) = \(vp\v a (x))\ 2 = \ipai3{.x)\ 2 from an initial 
state |^a(0)) to final state \v@) can then be brought into the following form 

.±fipaa\ 1 f-Am 2 cos29-2E u V NCB Am 2 sm29 \ (ijj aa \ 

% dx \ipab) ~ 4E U V Am 2 sin2# Am 2 cos 29 + 2E U V NCB ) \ijj ah ) ' W 

Here, E u is the neutrino energy and 9 = 9^4 and Am 2 = Am\ k with k — 1, 2, 3 for a = e, /x, r 
are the vacuum mixing angle and the mass squared difference between the new massive state 
z/ 4 and SM neutrinos ^1,2,3, respectively. The matter induced mixing angle 9m reads 

tan2# 
tan 29 M - { - ^-—^-^ ^ — , (10) 



and, since sign (Vjvob) — Qb, resonant flavor transitions for v h could be possible for q b = +1 
and 9 > 7r/4 or for qt> = —1 and 9 < n/A (and vice verse for F&.) We note, however, that the 
efficiency of a transition has a separate dependence on Am 2 , unrelated to Eq. fllOp . as the 
matter-induced oscillations seize to occur in the limit of Am 2 — > 0. From ([TO]) we find the 
necessary condition for which NCB effects are least likely to play a role, 

Am 2 cos 29 « l(T 4 eV 2 x (—??—) (— ) ( -^) ■ (11) 

VioMevy viooy Vg/ cm v 

In this work we focus on a parameter region which obeys this limit. A discussion for larger 
values of Am 2 is beyond the scope of this work. 

2.2 Solar v h flux 

Let us consider a scenario in which the baseline of z/& oscillation L osc is on the order of the 
earth-sun distance, L = 1AU ~ 1.5 x 10 8 km. This "just-so" choice of parameters suggests 
a canonical mass squared difference, 



°' 5X ^^ JllOMev'- (12) 



^OSC 



"Flavor" eigenstates v a i (a = e, /x, r, b) are found from mass-eigenstates uul by inversion 
of (jlj), v a L = J2k Uak^ki, and their evolution is obtained by solving 

i^=U% n = ^(UM 2 d U^ + A). (13) 

Here \I/ is the vector of amplitudes for the flavor states, \1/ = (ip^ip^,^,^), Ai 2 d = 
diag(m 2 , m|,m|,m|) is the diagonalized mass matrix and the entries of A = di&g(A C c + 
A N c, A NC , A NC , A NC b) are related to the induced matter potentials via A x = 2EV X . In gen- 
eral, the baryonic neutrino flux at the Earth location is found upon numerical integration 
of (TT3|) from the production point r$ of v e in the solar interior to earth at distance L with 
initial condition \l/(ro) = (1,0,0,0), This could be a complex problem when matter effects 
are involved, but fortunately not for the region of the parameter space we are interested in. 

With a few simplifying assumptions the appearance probability at earth can be obtained 
analytically [12j. We seek access to the high energy end of the neutrino spectrum, E v > 
10 MeV, because scatterings of Vf, will then more likely be picked up by a detector. The 
largest flux in combination with high endpoint energy comes from the neutrino emission in 
the decay of 8 B. With 4 He being the most tightly bound light nucleus, hep neutrinos have 
the highest endpoint in energy but come with a flux which is smaller by three orders of 
magnitude. The 8 B and hep respective fluxes and endpoint energies are given by |32j . 

$ 8B = (5.69+^^) x 10 6 cm" 2 s -1 , £ max ,8B = 16.36 MeV, 
$he P = (7.93 ± 0.155) x 10 3 cm' 2 s~\ E maxMp = 1.88 MeV. (14) 

In the solution to the solar neutrino problem the MSW effect dominates the flavor evolution 
of the highly energetic part of the neutrino spectrum and neutrinos exit the sun mainly as v 2 . 



Therefore, we assume a preferential mixing of the new neutrino to z/ 2 , # 24 7^ 0, and neglect 
other mixings of z/ 4 to the active states. Choosing q^ > 0, no resonant flavor transitions to 
V), inside the sun appear and with a ballpark of Am 2 suggested in (j!2p the standard solar 
MSW solution remains in place. 

With these assumptions and using a tri-bimaximal mixing ansatz for the active states 
the following v\> appearance probability at earth for the high energy part of the 8 B (and hep) 
fluxes has been obtained in [T2] . 



P eb {L,E u )~ sin 2 (29 b ) sin 2 



Am 2 b L(t) 



(15) 



It is assumed that mass mixings among active components are larger than mixings with 
v\, so that one can address the diagonalization of the neutrino mass matrix sequentially. 
In this procedure, Am 2 and 9 b denote the associated effective mass-squared difference and 
mixing angle between z/ 2 and Uf,, respectively. The true vacuum mass eigenstates are uj = 
cos 6^2 + sin 9bUb and i>u = — sin 6^ + cos 9 b Vh and a phase builds up between i/ 2 exit from 
the sun and propagation to the detector at distance L, 



L(t) = Lq < 1 — ecos 



2vr(t - t ) 



(16) 



where e = 0.0167 is the ellipticity of the earth's orbit; the perihelion is reached at to ~ 3 Jan. 
In what follows, it will be convenient to introduce an effective interaction parameter Af e s, 



Af 2 s =Af 2 sin 2 (29 b )/2. 
In the limit of rapid oscillations this implies P e bG 2 B —¥ Af 2 s G 2 F . 



(17) 



3 Direct Detection 



In this section we provide a detailed investigation of current and future direct detection 
experiments. From the scaling ([I]) we expect that elastic scattering off nuclei in direct 
detection experiments constitutes one of the most promising avenues in the search for a 
solar, long-baseline flux of v h particles. 

The spin-independent elastic recoil cross section on nuclei obtained from ([6]) essen- 
tially resembles the one from neutrino-nucleus coherent scattering [5] with the replacement 

G 2 F (N/2) 2 -)■ G\A 2 [E], 



dOr, 



d cos 9* 



§^ 2 a 



COS0*) 



'U 



Here A is the atomic number of the nucleus and 9* is the scattering angle in the CM frame. 
Equation (fl8l) can be rewritten in terms of a recoil cross section, 



da c 
dE 



H 



r^2 

2tt 



^m^dql) 



1- 



(E D 



E 2 



E 



R 



m N 



cos 0*) 



(19) 



where E min = ^jE R m N /2 is the minimum energy required to produce a recoiling nucleus 
of mass m ^ and kinetic energy Er. Here we have now included the nuclear form factor 
suppression F 2 (|q|) for scatterings with three-momentum transfer q, and in our numerical 
evaluations we use the Helm form factor parametrization [33] with the nuclear skin thickness 
of0.9fm. 

For the sake of comparison to the simplest case of spin- independent scattering of DM 
on nuclei, we can evaluate the total elastic scattering cross section of v\, on nuclei (at zero 
momentum transfer), 

ael = C 3l A 2 E 2 u ~ 1.7 x 10- 38 cm 2 x A 2 ( ^-) ' ( -r^j) ' . (20) 



Ti u \iooy v l0MeV , 

The coefficient in front of the second relation serves as a figure of merit when compared to 
the DM-nucleon cross section a n . Given that direct detection experiments have put upper 
limits on a n as low as 10~ 44 cm 2 [S] the coefficient in (|20p is sizable. However, a much more 
stringent cut off in Er makes it increasingly difficult for an essentially massless i/j, to scatter 
off heavier targets. 

The recoil spectrum arising from the solar flux of v\> will have to include an average over 
the neutrino energy spectrum dfi/dE u of neutrino source i, weighted by the Vb appearance 
probability and an overall flux modulation [Lq/ ' L{t)] 2 due to the earth's eccentric orbit, 



dR 

N 7 



dE R [L(t) 



dE v dE R 



Note that P e b depends on E v so that it has to be included into the average; $j is the integral 
flux given in (THj) and dfa/dE is normalized to unity, J dE dfc/dE = 1, and taken from [331135]. 
Nt denotes the number of target nuclei per unit detector mass, and in our computations we 
take the fractional isotopic abundances of each element under consideration into account. 
The rate exhibits a non-trivial time- dependence. The maximum of the overall flux is attained 
at the perihelion in early January. However, the integral in Eq. (T2~T|) constitutes an additional 
source of modulation which depends on the neutrino energy. It will have its strongest effect 
on the differential rate in the "just-so" regime of Eq. (IT2|) where L osc is on the order of the 
sun earth distance; we exploit this fact in the following section. 

The left panel of Fig. [T] shows the solar neutrino spectra of the various sources. In the 
right panel we compute the associated recoil spectra using (I2T|) for a germanium detector 
without threshold, with perfect energy resolution, and an exposure of 1 kgxyr. For simplicity 
and since it does not affect the argument we use P e & as given in (fl5l) for all fluxes, while 
strictly speaking Eq. (II 5p is only valid in the MSW regime. The deviation from MSW only 
affects the softest recoils, and can therefore be safely neglected in what follows. As can be 
seen, only 8 B and hep neutrinos reach out to values of Er > fewkeV where direct detection 
experiments become sensitive. Furthermore, 8 B neutrinos constitute the dominant part of 
the signal with hep giving a small correction only. 

The spectrum in (I2TI) is a theoretical one. To make contact with experiment we include 
effects from energy resolution, detector threshold and quenching of nuclear recoil energy. 
Details will be given when discussing the respective experiments. We start our discussion 
by considering the experiments with putative positive signal claims. 
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Figure 1: Left: Solar neutrino fluxes as a function of energy as taken from [3H |35, 36J. Right: 
Associated recoil spectrum in a (perfect) germanium detector with a total exposure of 1 kgxyr. 
Only 8 B and hep neutrinos reach out to values of Er where direct detection experiments become 
sensitive. 



3.1 DAMA 

The DAMA/Nal and its upgrade the DAMA/LIBRA experiment [37], situated in the north- 
ern hemisphere at the underground Gran Sasso National Laboratory (LNGS), were the first 
to report on a potential direct DM detection. The experiment uses large radiopure Nal(Tl) 
crystals to measure scintillation light resulting from nuclear recoils. Given that there is no 
other discriminating channel except requiring the candidate event to be a "single-hit", a 
relatively large overall count-rate of ~ 1 cpd/kg/keVee is observed. The presence of a posi- 
tive signal is inferred from the annual modulation of the residual count rate on the order of 
~ 0.02 cpd/kg/keVee over low-energy bins between 2 and 6 keVee once the average count rate 
per cycle is subtracted [131 E3] > f° r recent discussions on potential modulating backgrounds 

see [281132113011311132]. 

The modulation of the event rate has been observed over the course of more than a 
dozen annual cycles, collecting a cumulative exposure of 1.17tonxyr [131 113] • The null- 
hypothesis, i.e. a rate constant in time, has been excluded at the 8.9<r level. The residuals 
of the DAMA/LIBRA runs in consecutive bins between (2-4) keV, (2-5) keV and (2-6) keV 
are shown by the data points in Fig. [3J The DAMA signal is usually decomposed as 



S = S + S m cos [u(t — t )] 



(22) 



where Sq ~ 1 cpd/kg/keVee is the baseline rate of single hit events and S m is the modulation 
amplitude, 
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dK 



to 



dR 

dK 



(23) 



t +l/2yr. 



M 

M 

a 

o 

0) 

-a 

3 



o 

I 

o 



0.03 



0.025 



0.02 



0.015 



0.01 



0.005 



-0.005 



-0.01 



~ r ~ ' 'dama^oio h^T 



"ft 



■A4ff = 255 
Am 2 = 2.4x 10" 10 eV 2 




5 10 15 20 

E v [keVee] 

Figure 2: The data points show the DAMA modulation amplitude as reported in |14j in units of 
counts per day (cpd) per kg detector material and recoil energy. The solid line is the best fit from 
fft to the data. 



The measured phase t is reported to be compatible with the one expected from DM, t = 
152.5 days (June 2nd) with a period of one year, i.e. u = 2tt/(1 yr). The reported modulation 
amplitude is shown by the data points in Fig. [2j 

To see if Uj, provides a viable explanation of the DAMA data one can either fit the 
modulation spectrum (1231) or directly the time series of the residual rate. Considering the 
solar Uf, origin for DAMA, Eq. (12~T1) may not necessarily lead to a truly sinusoidal form 
of the signal as a function of time. In addition, t is n °t expected to be identical with 
the DM value of 152.5 days. At first sight, a direct fit of the time series seems therefore 
favorable. However, the reported residuals are binned in energy so that they only provide 
coarse-grained information on the recoil energy distribution. This, in contrast, calls for a fit 
of the modulation amplitude instead. We have implemented both approaches and discuss 
their results below. In addition, one can also attempt a joint fit of both data sets. This 
approach is complicated by the fact that the data sets are not independent. 

We start by fitting the modulation amplitude. Observable scatterings of Uj, occur on 
sodium only and no appreciable rate is expected for E v > 7keVee. The latter expectation is 
in accordance with what is seen in the data. Therefore, we only fit the first ten data points 
with E v < 7keVee in order not to bias the goodness-of-the-fit estimate. With the help of 
the usual y 2 function we obtain the following best fit values, 



DAMA S ri 



Am 2 b 



2.43 x KT 10 eV 2 , K 



cir 



255, 



Xmin/^d = 9.5/* 



(24) 



The minimum in y 2 is associated with a p- value of p = 0.3; n^ denotes the number of degrees 
of freedom. The result of this fit is shown by the solid line in Fig. El Confidence regions in 
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Figure 3: The data points show the DAMA/LIBRA reported residual event rate for various energy 
bins as a function of time. The red line is the residual event rate associated with the fit to the 
modulation amplitude in Fig. [2j The dotted line is a fit of the sinusoidal function A x cos[cj(t — to)] 
with u) corresponding to a period of one year and a phase to as expected from DM. As can be seen 
the Vb signal is approximately out of phase by one month. For a quantitative statement see main 
text. 



Ami an d A4ff are constructed by demanding, 



X 2 (Am 2 b ,Af e z) < xLn + AX 2 



(25) 



where Xmm * s given in (|24|) . We choose A% 2 = 9.21 which corresponds to generous 99% CL. re- 
gions. The choice results in the two disjoint gray shaded regions shown in Fig. [6j 

The above result looks promising. However, in contrast to the DM case one has to check 
how well the time dependence of the rate is met. The resulting signals in the various energy 
bins from the best fit values (1241) are shown by the solid (red) lines in Fig. [3j For completeness 
we also show by the dotted lines fits of the data by the sinusoidal function Ai cos [u(t — t )] 
with period 1 yr and to = June 2 nd as expected when the origin were due to DM scatterings. 
As can be seen by eye, the v\, signals seems to lag behind by approximately one month. 
Thus, the best fit corresponds to a phase inversion with a maximum rate at the aphelion 
with to ~ July 5th. For example, using the (2 — 4) keVee residuals, the best fit values of the 
modulation spectrum fl24l) . one finds x 2 / n d = 101/41 with p = 5 x 10~ 7 for the residuals. 
This points towards a very poor description of the full data. 

We can try to improve on the above situation by directly fitting the residual rate. This 
is an important check, since the time dependence of Eq. (I2T1) is non-trivial. Can we find a 
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corner in the considered parameter space in which one can alleviate the tension in the phase 
of DAMA and the i/& signal? From the (2 — 4) keVee data we obtain as best fit, 

DAMA residuals: Am 2 = 2.45 x l(r 10 eV 2 , Af cE = 183, xLJ n d = 73.2/41. (26) 

Though the fit fares slightly better on the residuals with p ~ 10~ 3 , this does not ameliorate 
the situation sufficiently. Moreover, the smaller value of J\f e s now somewhat under-predicts 
the modulation amplitude. Finally, even when we perform a joint fit (neglecting potential 
covariances of the data sets) we do not find any substantial improvement. We conclude, that 
even though the DAMA modulation amplitude is fit rather nicely, the time series of events 
speaks against the u b interpretation. 

As a final remark, we comment on the sodium quenching factor. For our analysis above we 
used Q = 0.3 in the conversion to electron equivalent recoil energy, S v (keVee) = QE R (keV). 
New measurements J43] seem to indicate 1) lower values Q ~ 0.15 and 2) a stronger energy 
dependence as previously thought. This has important implications for light DM as well as 
for the Ub hypothesis. We find that Q = 0.15 moves the DAMA regions in Fig. [6] towards 
larger values of J\f e s which are already excluded by all the other considered null searches. 
The situation then becomes more similar to the one already witnessed for DM. 

3.2 CoGeNT 

The CoGeNT experiment is a low-threshold nuclear recoil germanium detector situated in 
the Soudan Underground Laboratory. The latest data release is based on 442 live days taken 
with 0.33 kg target [16j [T7]. An unexplained exponential rise of the signal at lowest recoil 
energies 0.5-1 keVee is observed. The origin of it is unknown and has lead to the speculation 
that DM with a mass in the ~ 8 — 10 GeV range may be the cause of it. For spin-independent 
DM-nucleus scattering, the excess requires a cross section asi ~ 10~ 40 cm 2 . Such values are 
challenged by the null-result of XENON100 and by the low-threshold analysis of CDMS-II. 
However, more recently the collaboration has identified a source of surface-background events 
which may lead to a revision of the signal strength in the low-recoil bin 0.5-1 keVee. In the 
following we will investigate the possibility that the observed excess may be due to scattering 
of v h in the detector. We will also account for the possibility that the collaboration's results 
could be revised in the near future 



In addition to the signal-rise below 1 keVee the data also appears to be annually modu- 
lated in the 0.5-3.2 keVee bracket. The observed event rate peaks in mid-to-late April (2010) 
with a modulation amplitude of ~ 16%, most pronounced between 1.4-3.2 keVee [T7]. The 
latter behavior is neither expected from DM scatterings nor could it be explained by v\, 
scatterings since the recoil spectrum arising from 8 B neutrinos is cut off for E v > 1.4keVee. 
We also note that the modulation of nuclear recoil events in Ge in that energy regime has 
recently been challenged in a dedicated analysis by CDMS [44]. We will therefore not further 
address the potential modulation of the CoGeNT signal and await further data. 

Cosmogenically induced radioactive background has to be subtracted from the CoGeNT 
data before fitting the exponential excess. The radioimpurities in the crystal have been 
identified by the collaboration, with the most prominent ones given by the electron capture 
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Figure 4: Recoil spectrum from the 442 live-day run of the CoGeNT experiment. The black (gray) 
data points show the signal after (before) subtraction of the cosmogenic radioactive background. 
The solid line is a fit to the black data points. It decomposes into the contribution from Vf, (dashed 
line) and the contribution of a constant background (dotted line.) 



decays of 68 Ge and 65 Zn centered at 1.3keVee and l.lkeVee, respectively. From a fit of 
observed K-shell electron capture peaks seen in the high energy data and from the expected 
ratio of L-to K-shell decays one can subtract the low-energy L-shell background in the 0.5- 
3.2keVee window. We follow [45] in the subtraction and collect the time-stamped events in 
O.lkeVee bins. The result of this procedure can be seen in Fig. [4] as the difference between 
gray (with peaks) and black (peaks subtracted) data points. 

Nuclear recoil energies on germanium have to be converted into the measured ionization 
signal. We employ a Lindhard-type, energy dependent quenching factor, i? v (keVee) = Q x 
.E'ft(keV) 1 " 1204 with Q = 0.19935 [45] and account for a finite detector resolution by convolving 
the recoil signal with a Gaussian of width a 2 = (69.4 eV) 2 + 0.858 eV x ^(eV) [46J. Finally, 
the efficiency of the detector is provided in Fig. 1 of [17J. 

When fitting the CoGeNT excess at low energies we follow two approaches: In the first 
case we seek an explanation of the excess exclusively in terms of v b scatterings on Ge, allowing 
only for a constant background contribution. In the second case we relax the assumption 
on the background and allow, in addition, for an exponential background component, A x 
exp(— BE V ), with coefficients A and B to be determined in the fit. Clearly, in the latter 
approach a ^-induced contribution may not even be necessary as the excess resembles an 
exponential shape. This is therefore the most conservative way to treat the data in terms of 
new physics and it will show us the "compatibility" region in the (J\f e g,Amb) parameter space. 
As mentioned above, the CoGeNT excess is likely to be revised by the collaboration. In the 
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second approach, the additional exponential background is here to mimic that possibility 
without quantifying its (yet unknown) concrete strength. 

We first take the reported CoGeNT excess at face value and fit it by i/& scatterings on 
Ge together with a constant background rate. The best-fit values are 

CoGeNT: Am\ = 1.76 x l(r 10 eV 2 , A/" eff = 228, X*/n d = 33.6/24, (27) 

which corresponds to p ~ 0.1. The background rate is Co = 3.36cpd/kg/keVee. A finer- 
grained bin size improves the goodness-of-fit to x 2 / n d = 47.6/47 with p = 0.45. We consider 
this a very good description of the data. One should keep in mind that the subtraction of the 
cosmogenic background has uncertainties itself which are not accounted for in the errorbar. 
Figure H] shows the spectrum obtained from (|27|) . The dashed line shows the signal from i/& 
only and the solid line includes the constant background. Below E v < 0.5keVee the detector 
efficiency decreases rapidly, which explains the turn-off of the scattering signal. 

Figure |6] shows the inferred 99% C.L. regions in the (Am^Weff) parameter space as 
labeled. We use Ax 2 = 9.21, i.e. we treat the constant background rate as a nuisance 
parameter. This is equivalent of using a profile likelihood to infer the confidence regions. Two 
isolated islands are visible as labeled. The thin gray solid line labeled "CoGeNT hull" which 
touches the CoGeNT regions from above is obtained by allowing an additional exponential 
background in the fit (see discussion above.) Without further knowledge of the strength 
of this background, the full region below the line then becomes viable. Whenever the i/& 
signal becomes too weak, the background takes over in producing a viable fit. The general 
expectation is that once the collaboration revises their statements about the strength of the 
exponential rise, the CoGeNT favored regions will move vertically downwards, but at this 
point it is impossible to speculate by how much. In conclusion, we find that i/& can provide 
an excellent explanation to the CoGeNT data. 

3.3 CRESST-II 

The CRESST-II experiment [47], situated at LNGS, has recently presented their results from 
their DM "run32" with a total of 730 kgxdays effective exposure between 2009-2011 [21J. The 
analysis has been carried out using data collected by eight CaW04 crystals which measure 
heat and scintillation light resulting from nuclear recoils. The calorimetric phonon channel 
allows for a precise determination of the recoil energy with a resolution better than 1 keV. 
Nuclear recoils are again quenched in scintillation light. This is a virtue as it allows for a 
discrimination against e~ and 7 induced events. Moreover, the quenching factors of Ca, W, 
and O differ. To a limited degree, recoils against the respective elements can therefore be 
distinguished. 

The analysis [21] finds an intriguing accumulation of a total of 67 events in their overall 
acceptance region between 10-40 keV, shown by the solid line in Fig. |5j The low-energy 
threshold of each detector-module is determined by the overlap between e/7- and nuclear 
recoil band. Allowing for a leakage of one background e/7-event per module distributes the 
individual detector thresholds between 10.2-19 keV. Whereas e/7-events are a well control- 
lable background, the experiment suffers from a number of less well-determined sources of 
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spurious events. To assess how well the observed events can be explained in terms of new 
physics makes the modeling of such background — unfortunately — unavoidable. 

In the following we briefly mention each of the known background sources and outline 
our treatment of them: 

1. As alluded before, the thresholds of the detector modules are chosen such that a leakage 
of a total number of 8 e/7-induced events into the nuclear bands are expected. We find 
the energy distribution of these events by digitizing and binning the corresponding line 
from Fig. 11 of [21] • 

2. Degraded a-particles from radioactive contamination in the clamping system holding 
the crystals can be misidentified as nuclear recoils once their energy falls below 40 keV. 
A sideband analysis above that energy indicates that the distribution in recoil energy is 
flat. This allows to estimate the number of a-events in the acceptance region for each 
detector module and which is provided in Tab. 2 of [21]. We follow this prescription 
which yields a total of 9.2 events. 

3. Related to the previous source, 210 Pb a-decays from radioactive lead on the clamps 
holding the crystal and with the a-particle being absorbed by non-scintillating material 
constitutes another source of background. The peak at 103 keV full recoil energy 
in 206 Po is clearly visible and a fit of it allows to infer the overall exponential tail 
distribution in the acceptance region below 40 keV. The low-energy tail is due to 206 Po 
that is slowed down in the clamps before interacting in the crystal. We estimate the 
radioactive lead contamination of each detector module from the number of observed 
events in the reference region above 40 keV. This yields a nominal background of about 
17 events. 

4. Finally, low-energy neutron-nucleus scatterings in the crystals is a well known source 
of background. These neutrons can be produced by in-situ radioactive sources as well 
as by cosmogenic muons in- and out-side the detector housing. Since neutrons tend 
to scatter more than once, some information on the overall flux can be obtained from 
the amount of coincident events in different detectors. Such multiple scatters also 
tend to wash out the initial spectral information. The authors of [21] parameterize 
the neutron spectrum by a simple exponential dN n /dE = A x exp(— E/E^ec) where 
Edec — (23.54 ± 0.92) keV has been obtained from a neutron calibration run with an 
AmBe source. The best we can do is assuming a uniform neutron flux in all detector- 
modules. With A = 1 one gets a total of about 9 events. When fitting v^ to the data 
we leave A as a free parameter. We observe that A is never too large once the goodness 
of the fit becomes acceptable. The reason is that the ^-induced spectrum experiences 
are relatively sharp cutoff for energies in excess of ~ 25 keV. Thus, the high-energy 
part of the acceptance region has to be entirely explained by background for which 
A = 0(1) provides the best fit. v\> scatters mainly on oxygen and calcium, so that 
there is no need to further dissect the neutron background as the latter also scatters 
to 90% on O HO. 



In [2T] the various nuclear recoil bands have not been resolved. Therefore, we compute 
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Figure 5: CRESST-II recoil spectrum. The solid line is the histogram of reported events in the 
730 kg x days run summing to a total of 67 events. The gray shaded (stacked) histograms show 
the best fit contribution from v^ (darkest shading) and the modeled backgrounds as labeled and 
explained in the main text. The spiky dashed (red) solid line is the unbinned v\, signal. 
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rate predictions for CRESST summing up all events in Ca, 0, and W. The fractional expo- 
sures of the individual detector modules lie within ~ 20% of a uniform one with value 1/8. 
We use the accurate values as provided by (48]. We account for a finite Gaussian energy- 
resolution in the phonon channel with a — 1 keV. Since the number of observed events n, in 
each of the bins is small, we fit the data by minimizing the Poisson log-likelihood ratio 
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where the sum runs over all bins and i/j is the sum of background and signal contributions; 
the last term is absent when n^ = 0. Confidence regions are directly constructed from (T2"8"D . 

Figure |5] shows the recoil spectrum induced by i/& as the (magenta) continuous and falling 
line. Summing all background contributions to the uj, signal the best fit parameters read, 

CRESST-II: Am 2 b =3x l(T 10 eV 2 , A/" eff = 49, x*/n d = 27.7/27, (29) 

with a p- value p = 0.48 under the approximation that Xp i n (1281) follows a x 2 distribution with 
rid — 27. The amplitude of the neutron background is found to be A = 1.23. Fixing instead 
A = 1 yields the same parameters (1291) with negligible degradation in \ 2 ■ Discontinuous 
jumps in the count rate when going from lower to higher recoil energies are due to the onsets 
of the various detector modules with increasing energy thresholds. Also shown as a stacked 
histogram are the modeled sources of background as labeled. 

Figure |6] shows again the 99% confidence regions in the (Am 2 ,A/" e ff) parameter space. As 
can be seen, the favored region stretches across the plane and the trend for larger values of 
Ami n °t shown in the plot can be easily be extrapolated. The CRESST region is compatible 
with the one inferred from the DAMA modulation amplitude. Once the CoGeNT excess is 
revised (see previous section), it is very likely that the resulting best fit region will overlap 
with CRESST too. CRESST spans a rather wide region in parameter space. The light yield 
distribution of the candidate events as a function of E v is not published and has thus not 
been accounted for. With eventual better knowledge of this quantity, the region is expected 
to shrink in a joint fit. In the v^ scenario, most of the scatters are on oxygen which in turn 
yields most scintillation light among the CaW04 constituents. 

It is also likely that the new CRESST data constrains larger values of N" c s- This is 
especially true given that the detector model with lowest threshold only observed one count 
between 10-12 keV. In order to put a constraint we use what has been termed "binned 
Poisson" technique in [49] : for one bin, an average number of events v = u s + z^b g consisting 
of u s signal and v^ g background events is excluded at a level 1 — «bm if the probability to 
see as few as n Q ^ s observed events is cubm- Since n b s is a Poisson variable, ctbm is given the 
lower tail of the Poisson distribution, «bm = Yln=o ^ n exp(— v)jn\. When dealing with more 
than one bin the level of exclusion 1 — a is given by 

l-a = (l-a b in) 7Vbin , (30) 

where N hin is the number of bins; a is the probability to see as few events as observed in 
at least one of the bins. For placing a constraint from CRESST we only use the bins from 
10-25 keV as those are the ones for which i/j, can give a contribution. More conservative 
constraints are obtained when assuming that there is no background. 
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Figure 6: Summary plot of direct detection favored regions and constraints in the parameters Amj? 



and A/" e ff at 99% confidence. Favored regions: the broad light shaded gray band shows the CRESST- 
II region. The two darkest islands are the regions in which the CoGeNT excess is explained. In 
presence of an exponential background contamination (e.g. due to "surface events"), the region 
below the thin gray line labeled as "CoGeNT hull" becomes in principle viable (see main text for 
details.) The two medium gray shaded islands indicate the regions in which the DAMA modulation 
amplitude is fitted; these regions as well as any other parameter choices however exhibit a tension 
in timing when compared to the DAMA residuals. Constraints: A/" e ff values above the respective 
lines are excluded (or seriously challenged.) The top constraint is the one from XenonlOO, the two 
degenerate ones below are obtained with the CRESST-II data and CDMS-II low-threshold data. 
The two dotted lines at the bottom show the constraints arising from the XenonlO low-threshold 
analysis with two different assumptions on the ionization yield Q y (see main text for details.) 
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3.4 Null-searches 

3.4.1 CDMS-II low threshold analysis 

The CDMS-II collaboration has published a low threshold analysis from the Soudan site using 
eight Ge detectors with a raw exposure of 241 kg days [50]. At the expense of discriminating 
power of e~/7 against nuclear recoils, a threshold of 2 keV was reached. This is an interesting 
analysis because it uses the same target material as CoGeNT. Indeed, for the DM case the 
results indicate a serious conflict between the two experiments. Therefore, it is important 
to check if the CoGeNT explanation is also challenged in the baryonic neutrino scenario. 

An exponential-like signal rise towards threshold with a maximum event rate of ~ 
1 cpd/kg/keV has been observed; see Fig. 1 in [50J. "Zero-charge" events from electron 
recoils near the edge of the detector are expected to yield a major contribution to the ob- 
served count rate. Given that the estimation of this background involves extrapolation and 
hard-to-control systematic errors, we obtain the most conservative limits on A4a by not sub- 
tracting this background. This follows the approach taken by the collaboration itselfjj Only 
the first three data bins covering 2keV < Er < 3.5 keV are sensitive to v^ scattering. On 
those bins we perform a "binned Poisson" exclusion, similar to the one explained in Sec. 13.31 
We correct for efficiency according to Fig. 1 in [SU] and use a Gaussian detector resolution of 
0.2 keV. In Fig. Owe show the resulting constraint at 99%C.L. Remarkably, CDMS-II does 
not challenge the CoGeNT-favored regions. 

3.4.2 SIMPLE 

An interesting direct detection experiment in the current context is SIMPLE [52], operated 
in the Low Noise Underground Laboratory in southern France. It uses a dispersion of 
superheated liquid droplets made of Ci 2 ClF 5 with an total active mass of 0.2 kg. Only 
nuclear scattering induces phase transition which results in bubble nucleation. The fact that 
the experiment contains mainly light elements makes it susceptible to //^-scattering. Notably, 
fluorine with atomic number A = 19 has a target mass fraction of ~ 60%. 

We use the results from Phase I of Stage II with 14. lkg days exposure [52]. A total 
of 8 events were observed with an expected neutron background of 12. We include this 
background in the derivation of an upper limit of AW as a function of Ami- We model 
bubble nucleation and heat transfer efficiency following [52] for which we use an energy 
threshold of 8 keV. For better overview, we did not include the constraint in Fig. [6j It is 
superseded by the ones from CDMS-II and CRESST-II but more constraining than XenonlOO 
as will be discussed next. 



x It has been speculated [51] that the CoGeNT and CDMS-II recoil spectra are indeed similar after 
correcting for a potential energy miscalibration. Given that the status of the CoGeNT recoil spectrum is 
uncertain itself, we do not follow up on that discussion in this work. 



3.4.3 Xenon experiments 

In this work we consider the results from the XenonlO and its upgrade, the XenonlOO 
experiment at LNGS [53j IB]- Currently, the most stringent constraint for spin-inelastic 
scattering for DM masses in the 50 GeV ballpark is the one from the XenonlOO [8]. The 
experiment also provides strong limits in the light-DM mass region 0(10 GeV). In the latter 
context, the low-threshold analysis |54j from XenonlO is of particular interest. 

The experiments use a prompt scintillation light signal (SI) and a delayed one from 
ionization (S2) to detect the nature of the recoil event. Both signals are measured in units of 
photo-electrons (PEs). Nuclear recoil energies are obtained from the respective signals via 

SI : E R = -Lfifi, S2 : E R = ^. (31) 

Only a small fraction of the deposited recoil energy is emitted in form of scintillation light. 
The crucial quantity is the scintillation efficiency £ e e which it determines the discrimination 
threshold of the experiment. Therefore, for XenonlOO we will mostly be interested in SI 
since a low-threshold analysis has not yet been published. Conversely, for XenonlO we focus 
on S2. For SI we use the measurements of [55] with a conservative extrapolation of C c s to 
zero at 2keV nuclear recoil. L y = 2.2 (XenonlOO) is the light yield in PEs/keVee obtained 
from 7-calibration. S e and S n are quenching factors for scintillation light due to electron and 
nuclear recoils, respectively. They are experiment-specific and depend on the applied drift 
voltage; S e = 0.58 and S n = 0.95 for XenonlOO. Q y is the ionization yield per keV nuclear 
recoil and ( is the measured number of PEs produced per ionized electron, ( = 20 (24) PEs 
for XenonlOO (XenonlO); we will comment further on Q y below. 

The SI detector acceptance for XenonlOO is found from the lines presented in Fig. 2 
of [8] together with a low-energy threshold of 4PE which corresponds to 8.4 keV nuclear 
recoil energy. For S2, the trigger efficiency is effectively 100% in both experiments. We take 
the Poissonian nature on the expected number of scintillation photons/ionization electrons 
into account. For example, for S2 one computes 



f dR 

/ dE R -—xPoiss(n e \u e (E R )), (32) 



dR f ,„ dR 
dn f 



where n e = E R Q y is the number of ionized electrons; a PMT resolution of 0.5 PE can be 
neglected when converting n e to S2; S2 = n e (. We use a hard cut off at E m [ Q = 1.4 keV 
following [5"4"] . 

We first discuss the constraint from XenonlOO. Reference [8] presents the results from 
a 100 day run with a fiducial detector mass of 48 kg. Three events were observed in the 
acceptance region 8.4-44.6 keV. We use Yellin's maximum gap method [56J to set an upper 
limit. The resulting constraint is again shown in Fig. O Given that Xe is a rather heavy 
target and that the scintillation threshold is O(10 keV), the constraint on the parameter 
space is very week. Indeed, none of the favored regions is challenged. 

In contrast, a more stringent limit can be expected from the XenonlO low-threshold 
analysis in [51]. Discarding the scintillation signal allows one to lower the threshold to E R = 
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0(1 keV). After all cuts are applied and with a resulting effective exposure of ~ 6.15 kg days 
a mere number of seven events in the region 1.4 keV < Er < lOkeV are observed. The 
calibration of the nuclear energy recoil scale now solely depends on Q y [see Eq. (1HT|) ]. The 
problem in the analysis is that one requires Q y in a regime where no data is available and the 
extrapolation of that quantity is not well supported by theoretical expectations. Indeed, the 
adopted values in [51] by the XenonlO collaboration have been repeatedly disputed [571158] . 
Therefore, we adopt two different extrapolations of Q y . First we use the solid line in Fig. 2 
of [M] and second we employ more conservative choice discussed in [5B] in the estimation of 
the energy scale. 

The resulting constraints are again presented in Fig. In contrast to all other null 
searches, XenonlO severely challenges the entire region of interest with A4ff ^ 20. We 
caution the reader that it is also the constraint with the largest uncertainty. We will further 
illustrate the strong sensitivity of the results on Q y below when considering a prospective 
XenonlOO low-threshold analysis. 

3.5 Future sensitivity in Direct Detection 
3.5.1 XenonlOO low threshold 

In the previous section we have seen that a low-threshold analysis in XenonlO can yield very 
stringent constraints. Therefore, it is conceivable that the collected 100 live days of data 
with the XenonlOO detector may offer another sensitive test for this model. Here we present 
projections for a XenonlOO ionization-only (S2) study. 

Once the prompt scintillation signal SI is discarded, the goal lies in lowering the threshold 
in S2 as much as possible. In [8] the XenonlOO quoted software threshold is 300 PE which 
corresponds to 20 ionized electrons. This is a factor of five larger compared to the XenonlO 
low-threshold analysis. The gain in exposure by about an order of magnitude somewhat 
compensates for this since the increase in detector mass does not seriously affect the extrac- 
tion efficiency of ionized electrons. However, an air-leak during the run introduced unwanted 
Kr contamination at the (700 ± 100) ppt level. The associated rate in the electron recoil band 
reads, 

R Kv < 22 x 10~ 3 cpd/kg/keVee, (33) 

and is expected to be homogeneously distributed over recoil energy. Clearly, a putative v h 
signal must overcome this background. We remark in passing that for future runs (133]) will 
diminish since the 85 Kr concentration is being continuously reduced by cryogenic distillation. 



Figure [7J shows the sensitivity to z/& with the current dataset of 100 live days for Am 



2.5 x 10~ 10 eV 2 and A4fi = 100. We present the recoil spectrum as a function of actual 
detected S2 in units of PE. The solid lines show the 8 B and hep neutrino spectrum as 
labeled. The vertical dashed line indicates the XenonlOO threshold and the horizontal one 
is the rate (I3"3"|) from Kr contamination. As can be seen, a very large signal is possible. The 
dotted line is again a 8 B spectrum but this time with a calibration scale following [58J. As 
a consequence the signal falls entirely below the threshold. Even a small change AEr in 

20 



10000 



1000 



XENON100, 100.9 live days 




400 500 600 

S2 (PE) 



Figure 7: Projection for a XenonlOO low-threshold analysis for an exemplary parameter choice 
Amj = 2.5 x 1CP 10 eV and N e g = 100. The x-axis gives the ionization signal S2 in units of PEs. 
The horizontal and vertical dashed lines show the maximum rate from radioactive Kr decay and 
the S2 software threshold of the detector, respectively. The solid lines are the tvsignals from 8 B 
and hep neutrinos as labeled. The vertical arrow at 700 PEs indicates the current threshold of the 
SI scintillation signal. The dotted line shows again 8 B neutrinos for a calibration scale following 
[58] instead of [54] . This highlights the severe sensitivity on the extrapolation of Q y . 
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Figure 8: Predictions for a COUPP 60 kg bubble chamber for Ara^ = 3 x 1CT 10 eV and N e g = 60 
as a function of bubble nucleation threshold with one year of exposure. The solid lines are from 
top to bottom for the 8 B and hep fluxes of v\,. For comparison, a signal from a 10 GeV DM particle 



with spin-independent nucleon cross section of a n = 10 cm is shown. 



the nuclear recoil calibration has a large effect since AS2 ~ £AEr. This illustrates 1) that 
the XenonlO constraints in the previous section should be viewed with care and 2) that 
without further experimental insight into Q y in ( 13TJ) a conclusive prediction for XenonlOO is 
not feasible. 

3.5.2 COUPP 



As already mentioned in Sec. l3.4.2| detectors which employ fluorocarbon compounds as target 
material are attractive because of their favorable kinematics in contrast to heavier targets. 
A large scale experiment of this type is the COUPP 60 kg bubble chamber currently in the 
progress of moving into the SNOLAB underground facility [59]. It uses a superheated CF 3 I 
liquid with temperature and pressure adjusted such that only nuclear recoils set off bubble 
nucleation events. It is a counting experiment for events above adjustable threshold without 
a priori insight into the recoil energy distribution. We assume an exposure of 1 yr together 
with a detector efficiency e = 0.7. 

Light target nuclei make COUPP particularly attractive for the searches of light WIMPs 
and v h . In Fig. [8] the integral signal of u b for Am 2 = 3 x 10~ 10 eV 2 and A/" e fr = 60 as a function 
of detector threshold is shown. The two solid curves from top to bottom correspond to 8 B and 
hep neutrinos, respectively. Already below a threshold energy of ~ 20 keV the 8 B flux induces 
a clear signal. Limited insight into the energy distribution should be possible by varying 
the rather "steplike" detector threshold. In particular, with a multi-year exposure (or larger 
values of A r cr r) the crossover from the 8 B to the hep neutrino spectrum may be observable. 
More importantly, the variation of the threshold provides discriminating power between a 
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putative DM signal and v\>. The dashed line in Fig. [8] shows the integral event rate for a 
lOGeV DM particle with spin-independent WIMP-nucleon cross section of a n = 10 -41 cm 2 . 

4 Neutrino searches 

Here we are going to consider the elastic scattering of v\, from the sun in solar neutrino 
experiments. The NC channels in those experiments are not necessarily sensitive to this 
class of new physics, given that the associated inelastic reactions exhibit the scaling ([Q). As 
shown in [T2], the NCB interaction does not yield an observable rate for the D-breakup at 
SNO. There is also the possibility of inelastic excitation of 12 C with subsequent emission of 
a 4.44 MeV 7. The analysis of inelastic processes falls outside of the scope of the present 
study. 

The experiments which are capable of detecting the elastic NCB signal employ hydro- 
carbon scintillators, namely, Borexino and KamLAND as well as the upcoming experi- 
ment SNO+. The dominant background at lowest energies comes from 14 C contamina- 
tion of the mineral oil which decays with a maximum /3 energy of Q = 156 keV. Other, 
less prominent backgrounds by long-lived isotopes at lowest energies are the /3-decays of 
85 Kr (Q = 687keV) and 210 Bi (Q = 1.16 MeV). Unfortunately, the decays from 14 C pre- 
vent sensitivity to Uf, in current detectors. For example, the Borexino experiment which 
uses dedicated mineral oil with low residual 14 C content measured its concentration to be 
14 C/ 12 C ~ 2 x 10 -18 . Though this is a seemingly small ratio, it translates into a rate of 
approximately r ~ 6 x 10 4 events/day /ton below 0.2 MeV. Therefore, in the following we 
restrict our analysis to a future possibility when the level of 14 C is much reduced so that one 
may hope to gain sensitivity to Uj,. 

The scattering of 8 B neutrinos on protons can produce a recoil energy of Er < 0.5 MeV. 
However, the proton recoil is quenched and the v\, signal may not show itself above the 
14 C peak. In organic scintillators Birk's law provides a phenomenological description of the 
scintillation light yield per unit path length 



dL dE/dx 

dx 1 + ksdE/dx 1 

where k B is Birk's constant and dE/dx is the ion stopping power in the material. The 
formula interpolates between the limiting cases of low energy losses (no quenching), dE/dx <C 
kg 1 with linear dependence of the light output dL/dx ~ L^dE/dx and high energy losses 
(quenching), dE/dx ^> kg 1 for which saturation occurs, dL/dx ~ L /k B . Equation (13"4"|) 
suggests the (non-linear) relation between recoil and quenched energy, 

[ Er dE 

v ~ J l + k B dE/dx' ^' 

In the following we use Borexino, representative for other liquid scintillator experiments. 
It has a a fiducial detector mass of 0.278 kton filled with the scintillator pseudocumine, 
C 9 H 12 , with a mass density of p = 0.88 g/cm 3 . We employ the SRIM computer pack- 
age to obtain the stopping power dE/dx of protons in pseudocumine. Quenched energies 
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Figure 9: Liquid scintillator neutrino detector filled with pseudocumine, C9H12, and modeled after 
Borexino. Shown are the signals for v\, for J\f e g = 100 and Am 2 = 2.5 x 10~ 10 eV . The spectrum 
is largely dominated by background decays of 14 C for which a contamination in carbon of one part 
in 10 -23 has been assumed. Solid (dotted) curves are with (without) detector-resolution. Proton 
recoil energies are quenched; the quenching of electrons at lowest recoil energies has been neglected. 



are then obtained from ( 1351) with ks = 0.01 cm/MeV. The scintillation light yield is ap- 
proximately 500 PE/MeV and, for simplicity, we assume a Gaussian energy resolution with 
a = 0.045 MeV y/E^ [61] where E v is in units of MeV. 

We simulate the 14 C background spectrum as follows. The decay rate from 14 C decay in 
C9H12 is given by 



dT 



ii 



#14 9N A df u 3.1 x 10 4 



dE ti/2 In 2 Mp C dE day x ton \ 10 
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11 



1-18 
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dE 



(36) 



Here, R 14 is the ratio of 14 C/ 12 C, t 1/2 = 5730 yr is the half-life of 14 C, N A is Avogadro's 
number, and Mpc = 120.2 g/mol is the molecular weight of pseudocumine. The /3-spectrum 
is given by 



df 



11 



dE N 



- x Pe E e (E - E e ) 2 F(Z, E e )C(E e 



(37) 



where p e and E e are the electron momentum and total energy, Eq is the total endpoint 
energy and F(Z, E e ) is the Fermi function for 14 N; for the shape factor we use C = 1 — 
0.7i? e /MeV [53]. A^ is chosen such that the spectrum is normalized to unity. Electron re- 
coils are somewhat quenched by roughly (10 — 40)% for (100 — 10) keV [63J. For simplicity 
we neglect this complication here as well as more serious resolution effects near the detector 
threshold. In a realistic scenario, the latter can affect the spectrum below 100 keV substan- 
tially, see e.g. [63j. Here we are merely concerned with the question of how low the 14 C 
content needs to be in order to gain sensitivity to u h recoils. 
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In Fig. [9] we compare the 14 C background to a solar i/& signal with parameters A4ff — 100 
and Am 2 = 2.5 x 10~ 10 eV 2 . The contamination has been fixed to R\^ = 10 -23 which is 
already 5 orders of magnitude below the level for the Borexino detector, i? 14 (Borexino) ~ 
2 x 10~ 18 [6l]. There is little hope that a signal would show its presence above the 14 C- 
decay endpoint energy of 156 keV. The detector resolution has the effect of smearing the 
cutoff out to larger values of energy, hence burying the neutrino signal. This can be seen 
by the difference between the dotted and solid lines. At lowest recoil energies E R < 60 keV 
the 8 B signal dominates, but this may be a region which will not be explorable in large 
mass scintillator experiments due to their detector thresholds. Perhaps the only reasonable 
hope for using carbon-based scintillators is the isotopic purification of not excessively large, 
0(100 kg), quantities of carbohydrates with the eventual setup similar to the prototype of 
the Borexino detector 



5 Conclusions 

We have performed an extensive analysis of the model of v^ neutrinos that are sourced by 
the sun, and elastically scatter on nuclei in underground dark matter experiments. The goal 
of this study was to assess the viability of this model as the explanation for the reported 
anomalies, that are often interpreted as a possible dark matter scattering signal. Our findings 
are summarized in the master plot of Figure 6, that is a direct analogue of the WIMP-mass 
vs. scattering cross section plot for DM scattering. We are now able to conclude the following: 

• On the positive side for the model, the v^ scenario with an oscillation length compara- 
ble to the earth-sun distance and with an effective enhancement of the NCB current 
by O(100) is currently not seriously challenged by any of the existing experiments. 
On the contrary, this corner of parameter space provides a natural fit to the CoGeNT 
excess, and to the CRESST anomaly. Given enough uncertainty in the current status 
of the CoGeNT excess, and in the background-contaminated CRESST events, it is not 
difficult to see that the same regions of the parameter space of v\> model can be respon- 
sible for these anomalies. The null results of many experiments that challenge WIMP 
explanations of these anomalies (such as CDMS and XenonlOO), do not challenge z/& 
model. 

• On the negative side for the model, the values M c s — 200 which are necessary to fit the 
DAMA modulation amplitude are challenged by the CDMS-II low threshold analysis 
and by recent CRESST-II results. It should be noted however, that in the presence 
of background [HS], smaller values for A4ff become viable again. This requires larger 
modulation amplitudes of the signal, but this can presumably be achieved. Also the 
phase of the predicted modulation signal is deviant from the DAMA data. Even though 
one can achieve a phase reversal and have a maximum in the v\, scattering rate in early 
July, this is still far away from the DAMA phase, which corresponds to a modulation 
maximum in late May-early June. We do not know how to "correct" the model per 
se for this residual discrepancy. Further tension for the model may arise from a large 
number of events predicted for the "ionization-only" signal at XenonlO. At face value, 
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the parameter space of the model can potentially be constrained down to M ~ 10, but 
the severity of constraint may well be mitigated by a poorly known energy calibration 
at those lowest recoils. 

• Future prospects for probing ^-scattering look reasonably bright. In contrast to many 
DM models, the v\> scattering pattern is fixed and we can make definite predictions 
as functions of only two parameters. For example, these predictions show that the 
COUPP experiment may be particularly sensitive to the i/& scattering signal. One can 
also conclude that a re-designed iteration of the CRESST-II experiment with reduced 
backgrounds will likely be very sensitive to the v^ model. Furthermore, a reduced 
electron-like background in the XenonlOO experiment will be able to probe deep into 
the parameter space once the uncertain values for Q y and/or £ eff are clarified. 

Finally, as this paper was readied for the submission, a new preprint appeared [66J that 
examines a similar set of ideas. It expands the set of interesting mediation mechanisms 
beyond NCB to e.g. B — L and "massive photon" forces with extra-light mediators in the 
sub-MeV mass range. Unlike the Vb model, such modifications may already be under strong 
tension from astrophysical and cosmological constraints. We plan to return to //^-related 
signatures in astrophysical and cosmological settings in future work. 
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